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ABSTRACT 



o 

u 

O i Mixing of blackbodies with different temperatures creates a spectral distortion which, at lowest 

q order, is a y-type distortion, indistinguishable from the thermal y-type distortion produced by 

+3 ■ the scattering of CMB photons by hot electrons residing in clusters of galaxies. This process 

occurs in the radiation-pressure dominated early Universe, when the primordial perturbations ex- 
cite standing sound waves on entering the sound horizon. Photons from different phases of the 
CM sound waves, having different temperatures, diffuse through the electron-baryon plasma and mix 

together. This diffusion, with the length defined by Thomson scattering, dissipates sound waves 
' and creates spectral distortions in the CMB. Of the total dissipated energy, 2/3 raises the average 

00 ' 

^v^j ' temperature of the blackbody part of spectrum, while 1/3 creates a distortion of y-type. It is well 

ly-^ ' known that at redshifts 10 5 < z < 2 x 10 6 , comptonization rapidly transforms y-distortions into 

■ a Bose-Einstein spectrum. The chemical potential of the Bose-Einstein spectrum is again 1 /3 the 

value we would get if all the dissipated energy was injected into a blackbody spectrum but no 
extra photons were added. We study the mixing of blackbody spectra, emphasizing the thermo- 
dynamic point of view, and identifying spectral distortions with entropy creation. This allows us 
to obtain the main results connected with the dissipation of sound waves in the early Universe in 



a very simple way. We also show that mixing of blackbodies in general, and dissipation of sound 
waves in particular, leads to creation of entropy. 

Key words, cosmic background radiation — cosmology:theory — early universe — 

1. Introduction 

The Cosmic microwave backgr ound (CMB) has a s pectrum which is close to a blackbody with very 



high precision. COBE/FIRAS (iFixsen et alJll996l) constrained possible departures from a perfect 
blackbody to be < 10~ 5 for y-type distortions, and the chemical potential was limited to |//| < 
9 x 10~ 5 . This observation immediately has an important consequence: there was a time in the 
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history of the Universe when matter and radiation were in complete thermodynamic equilibrium 
with each other and any energy release at z < 2 x 10 6 was small. 

In addition to being an almost perfect blackbody, the CMB is also isotropic, with anisotropics 
of less than 10~ 3 . The dominant component of the anisotropy is the dipole caused by our pecu- 
liar motion. The anisotropies other than the dipole are at the level of 10~ 5 and originate from 
the primordial fluctuations imprinted in the initial conditions of the present expanding Universe. 
The radiation field seen by any observer (or electrons in the hot gas in a cluster of galaxies/early 
Universe) will thus consist of blackbodies having temperature that differs as a function of observa- 
tion direction. Scattering by electrons, averaging of the smaller scales within the beam of the tele- 
scope or explicit averaging by a cosmologist, will thus inevitably mix these b lackbodies together. 



Mixing of blackbodies was first studied in detail by (IZeldovich et al 



1972) who showed that it 



creates a .y-type distortion, indistinguisha ble from the y-type distortion c reated by interaction of the 
blackbody photons with hotter electrons (IZeldovich & Sunvaevlll969l) . IChluba & Sunvaevl ( 120041) 



proved that at lowest order, this result is valid for arbitrary temperature distributions and not just 
the Gaussian ones. Comptoniza tion of this initial y-type distortion converts it into a Bose-Einstein 
spectrum or a /U-type distortion ( Zeldovich & Sunvaevl 1969[lchluba et al.ll2012bl) fory > 1, where 
the Compton y-parameter is defined as the following integral over time t: y — J« e cr T c^fdf, n e 
is the electron number density, cr T is the Thomson cross section, c is the speed of li ght, kg is the 



Boltzmann's constant, T e is the el ectron temperature, and m e is the mass of electron dKompaneets 
Zeldovich & Sunvae vl 9l urpose of this paper is to study the mixing of blackbodies 



1956; 



from a statistical physics point of view. Amazingly, using basic thermodynamic relations, we arrive 
at the main results connected with the dissipation of sound waves in the early Universe due to shear 
vi scosity and thermal co nduction in a very simple manner. Some of these aspects are also discussed 



Chluba et al 



(2012b). 



2. Entropy of a Bose gas 

The entropy density of a general distribution of a Bose gas (photons), which may or may not be in 
equilibrium, is given bv Landau & Lifshitzll98ol) 



" 87r (^f) J" x2[(1 +n)ln(1 +")-" ln (")] dx (!) 

where = 8?r ^~ dv is the density of available states in the frequency interval dv, and for photons 
the degeneracy g = 2, h is Planck's constant, n = is the photon occupation number and E v 
is the energy density of photons per unit frequency. We have changed variables to dimensionless 
frequency x = hv/ik^T) defined with respect to a reference temperature T in the second line. For 
a blackbody spectrum at temperature T with n = rip\(x) = l/(e l - 1), we obtain S p \ = jrrm = 
|flRT 3 /A:B, where the last line defines the radiation constant «r. 

Now let us add a small spectral distortion 5n to the Planck spectrum so that the total occupation 
number is n — «pi( 1 + 6n/np{), where 5n /rip\ « 1 . Expanding the expression for entropy up to first 
order in 611 / «pi and ignoring the higher order terms, we get S — S p i + SS , where SS is the entropy 



[(1 + ri) ln(l + n) - n ln(n)] 
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density added/subtracted due to the spectral distortion and is given by, 



( 1 + npi) - ln(«pi)] 5n dx 



(2) 



where 5Q is the energy density in the distortion. Thus for small distortions, the classic thermody- 
namical formula for the change in the equilibrium entropy due to addition/subtraction of energy 
remains valid even in non-equilibrium. We note in particular that, at first order in distortions, the 
entropy does not depend on the shape of the distortion but only on the total energy density in the dis- 
tortion. This means that the process of comptonization, which converts an initial y-type distortion to 
a /i-type distortion, in the absence of any additional heating/cooling, does not change the entropy 
of the photon gas at first order. The factor of ks in Eq. (|2]i just converts temperature from Kelvin to 
energy units to make entropy dimensionless. As an example, we can calculate the entropy produced 
during cosmolo gical recombination, when ~ 5 energetic recombination photons are produced per 



hydrogen atom (|Chluba & Sunvaev 



(IRubino-Martin et al 



2006). Adding up the energy in the recombination spectrum 

20081) . We get <J5 recombination ~ ^Grecombination/C^B^CMB) ~ 9 X 10~ 6 /cm 3 to- 



day, compared with entropy density of 1478 /cm 3 in the blackbody part of CMB with temperature 
Tcmb = 2.725 K. 

In the rest of the paper we will measure frequency and temperature in energy units and set 
c = h = &b = 1 unless explicitly specified in definitions of other constants. 



3. Mixing of blackbodies 

Since the blackbody radiation is an equilibrium distribution for photons it is described by a single 
parameter, the temperature T. In addition, if we specify a second thermodynamic quantity such as 
volume V, we can calculate any other thermodynamic quantity such as entropy or internal energy. 
The spectrum is described by the well known Planck function, and the energy density per unit 
frequency and total energy density integrated over frequency can be written as, 

E v = 8ttv 3 



■t - 1 



E = a s T*. (3) 



The number density of photons is given by N = b^T 3 where b& = l6nl ^J^ 3 \ f is the Riemann 
zeta function with f(3) ~ 1.20206. We can also calculate the entropy d ensity and it is given by 



S = 4E/(3T) = {Aj3)a R T\ For the CMB with Tcmb = 2.725 K+lmK dFixsen & MatheJI 20021) 
we have energy density Eqmb = 0.26 eV/cm 3 , number density A^cmb = 411/cm 3 and entropy 
density Sqmb = 147 8 /cm 3 . For simplicity, we will consider mixing of two blackbodies below, 
but the derivation and the results are trivially generalized to an ensemble of arbitrary number of 
blackbodies by just replacing the average of quantities over two blackbodies with the appropriate 
average over the whole ensemble. 
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3.1. y-type distortion 

If we mix blackbody spectra with different temperature, the re sultant spectrum is not blackbody and 



at lowest order the distortion is given by a y-type spectrum dZeldovich et al 



3B 



1972)H This can be 



seen immediately by Taylor expanding the photon intensity or equivalently the occupation number 
n of a blackbody at a temperature T + AT about the average temperature T and take the average 
(ensemble or spatial) keeping terms up to second order in AT/T, 

1 



(n Pl (T + AT)) = 
[nn(T) ■ i 1 + 
=n Pl (T) 



AT 



dn Pi 



1 



KfM(f) 



dln|T] 



In 



1 + 



dnpi 1 
T— + - 

dT 2 



AT 
~T~ 

ATY 



d 2 n P i 
d(MT]) 2 i 

T 4 d I dnn 
dT T 2 dT 



="pi 



1 + 



ATY 



+ -Y(x) 
2 V ' 



ATY 



(4) 



where we have used (^r) = 0. The first term in the last line is simply a blackbody with temperature 
Fnew = T[l + (AT/T) 2 ] > T, and 

e x + 1 



Y(x) = 



(e x - l) 2 \ e x - 1 



-4 



is the y-type spectrum with x - v/T dZeldovich & Sunvaev|[l969l) with the magnitude of distortion 
given by yy — \ ^(t") ^- The average spectrum with y-type distortion for average of two black- 
bodies with temperature T + AT is shown in Fig.Q] Figure [2] shows the same spectra but in terms 
of the effective temperature defined by n(v) = l/(e v/Tca( ^ - 1), which for small distortions can be 
written, at linear order, as T e g(x) = T + T [(1 - e~ x )/x] 5n/n. We note that in the Rayleigh-Jeans 
limit, intensity is proportional to temperature, and thus the average of blackbody spectra just gives 
the spectrum at average temperature T, and T e g — > T. An important property of y-type distortion 
is that it represents pure redistribution of photons in the spectrum due to addition of energy but 
conserves photon number. Thus J x 2 Y(x) dx — and the above mentioned decomposition of the 
spectrum into a blackbody part and y-distortion part is unique and independent of gauge/reference 
frame, at order Sn/n, in the following sense: The constraint J x 2 Sn(x) Ax — on the spectral dis- 
tortion 5n fixes the reference temperature of the blackbody part of the spectrum used to define the 
variable x and 5n as a function of the variable x is gauge independent^ 

Without loss of generality, let us consider the superposition of two blackbody spectra with 
temperatures T\ — T + AT and To — T - AT, with average temperature T. The average initial 
energy density, number density and entropy density of the two blackbodies is 



(5) 



^initial - r^R^l + ^l) 



^initial = + T l) 



a R T 4 
b R T 3 







1+6 










^Ar\ 2 


1+3 









>aaT 



> bpT 



1 4flR , , 4flR , 

"J initial — ^ ~~ ~ ^1 + 1 2 ' ~ ~ ~ 1 



1+3 



(f) 



4«r -> 
> —T 3 
3 



(6) 



1 We discuss different processes which can lead to the mixing of blackbodies in the following sections. 

2 In fact 8n(x) is also invariant to change in the reference temperature happening in the mixing of black- 
bodies, Sn(x) = 5n(x ncw ) + 0{[6n/n] 2 ), where x ncw = v/r ncw . 
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Fig. 1. Average of blackbodies of temperature T + Ar and T - AT creates a new spectrum marked 
'Average(Y)' which is different from blackbody spectrum marked with temperature T. The aver- 
age spectrum is just the usual y- type distortion with respect to a new blackbody at temperature 
T[l + (AT/T) 2 ], with the two crossing at x = 3.83. At redshifts z > 10 5 the average spectrum 
will comptonize to Bose-Einstein spectrum marked fi above. All three spectra, T[l + (AT/T) 2 ], 
Average(Y), and fi, have the same number density of photons. Average/y-type spectrum and Bose- 
Einstein//i-type spectrum also have the same energy density which is greater than the energy density 
in the blackbody spectrum with temperature T[\ + (AT/T) 2 ] by 1/3 of the initial energy density 
excess over that of blackbody with temperature T. We have used linear order formulae to calculate 
the y and fi distortions in the figure but used a large value of AT/T to make differences between 
different curves visible. 

We can calculate the final temperature of a blackbody having the same number density of photons 
as the initial average A^Mai- 



nal 



/ ^initial \ 

\ b R ) 



1/3 



1 



f) 



(7) 



This is exactly the temperature, r new , of the blackbody we got by averaging the intensity in Eq. ©. 
Thus all the initial photons go into creating a blackbody with a higher temperature. The entropy 
density of this new blackbody is also identical to the initial average entropy density S initial because 
number density and entropy density have the same T 3 temperature dependence. The energy density 
of the new blackbody is however given by, 



£nnal - CIrT 



final 



i 



<E 



initial ■ 



(8) 
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Fig. 2. Same as Fig.Q]but effective temperature T e ff (v) defined by «(v) = 1 / (e v ' Tes(y ' — 1) is plotted. 
We have used linear order formulae to calculate the y and fi distortions in the figure but used a large 
value of AT /T = 0.4 to make differences between different curves visible. Average temperature 
T = 2.725 K. For Bose-Einstein spectrum we thus have, at linear order, T e g(x) = Tbe - Tfi/x. For 
a general distortion Sn/n (including y-type ) we have T e ff(x) = T + T [(1 - e~ x )/x] bnjn. 



In fact we find 



£final - a R T 4 - 2/3 ^initial ~ CIrT 4 ^ . 



(9) 



This result can also be obtained directly by multiplying Eq. © by v 3 and integrating over frequency. 
Equation (|4]i also shows that the rest of the initial energy, equal to 

(AT\^ 1 I AT\^ c 

— I a R T 4 K- I — I J dxx 3 Y(x), (10) 

goes to the y-distortion wi th the magnit ude of the distortion given by yy = \ [^f \ ■ These results 
were recently obtained in 



(Zeldovich & Sunvaev 



Chluba et al 



(12012bl) using the Boltzmann equation. It is well known 



1969) that v-distortion of magnitude yy decreases the brightness tempera- 



ture of radiation in the Rayleigh- Jeans part of the spectrum by an amount -2yy which is equal to 
(4£) for the mixing of blackbodies considered here. But we have also increased the magnitude 
of the brightness temperature by same amount in the blackbody part of the spectrum, resulting in 
the brightness temperature in the Rayleigh- Jeans part for the total spectrum which is equal to the 
average temperature T of the initial blackbodies, as shown in Figs. Q] and [2] We can calculate the 
additional entropy in the final spectrum resulting from the mixing of blackbodies using Eq. d2]i, 

■ 2 



AS 



AE 



2a R T J 



AT 



(11) 



where AE =1/3 feinitial ~ grT 4 ^ is the energy in spectral distortion, i.e. deviation from the black- 
body spectrum. 
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To summarize, mixing/averaging of blackbodies leads to a new spectrum which has a black- 
body part at a temperature which is higher than the average temperature of initial blackbodies by 
Ar ne „ = T(AT/T) 2 . This new blackbody has the same entropy as initial average entropy and the 
same number of photons but less energy. The remaining energy density appears as a y-type distor- 
tion, which is just a redistribution of the photons of the new blackbody spectrum. The y-distortion 
part can be identified with the increase in the entropy of the system (Eq. (fTTTi), which is expected 
as there is an increase in disorder of the system. 



3.2. /u-type distortion 



We saw in the previous section that it is impossible to create a blackbody spectrum by mix- 
ing/averaging blackbodies of different temperatures. The reason is that there is too much energy 
density for the given number density of photons. However we can still create a Bose-Einstein spec- 
trum with a non zero chemical potential /j., which is the spectrum we would get if the photons could 
again reach equilibrium while conserving energy and number. It is straightforward to calculate the 
temperature and chemical potential of the Bose-Einstein spectrum, «be = 1 /(e v ' TBB+ft - 1), by equat- 
ing the initial average photon number A^iai and energy density initial to the photon number A^be 
and energy density £be in a Bose-Einstein spectrum, as is done in case of heating of CMB without 
adding additional photons ( Illarionov & Sunvaev 19751) . In the limit of sm all chemical potential, 
li <s 1 , with Tee = r new (l+?BE) ~ T (l + ?be + [Ar/T] 2 ) , ?be << 1 we have ( Illarionov & Sunvaev 



1975 



a 



^£^7^(1-1.1106//) 



a R T 



^initial ~ CIrT 



(ATV 



+ 4| 

2 



1+6 



#be*W e (1 - 1.3684^) 



= N ir 



b R T 



b R r 



"f 



1.1106a/ 



( /Ar\ 2 

l+3f BE + 3^— j - 1.3684a/ 



(12) 



The second equation, A^be = Mnitiah directly gives us the relation ?be = 1.3684///3 = yu/2.19. 
Solving the system of equations for ?be and ji in terms of AT/T, we get 



H = 1.4 



2 (^initial - flR^ 4 ) - 2.8 



tBE 



JL- 

2.19 



(AT 
ii — 



(13) 



This is exactly the result we will get if we add energy AE = | (initial - crT 4 ^ to a blackbody 
with temperature r ne „ while conserving the photon number. We can thus write the deviation of the 



3 Using 6<T(3)// 3 ~ 1.1106 and n 2 /(3I 2 ) ~ 1.3684, where f is the Riemann zeta function with £(3) ~ 
1.20206, and h = f x 3 n PI (x) dx = tt 4 /\5, I 2 = f x 3 n PI (x) dx = 2f(3). 
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yu-type spectrum from the blackbody with temperature r new : 

»be - »pi(r ne w) . . e x 

™ — r =ObeX-//)— — 7- (14) 

"PlU new) e x - 1 

The yu-type spectrum resulting from the average of two blackbodies is also shown in Figs. Q] 
and |2 An important point to note here is that the Bose-Einstein spectrum is uniquely fixed by 
energy density and number density constraints. In particular, the value of the chemical potential 
is independent of any reference temperature we may choose to define the dimensionless variable 
x. This is, of course, the same value of /u we would get if we just comptonize the y-type distortion 
of the previous section. The Bose-Einstein spectrum and the chemical potential has however more 
fundamental origins in statistical physics, compared to the y-type distortion, the shape of which 
originates in the Compton scattering process or sum of Planck spectra. The //-type results of this 
section can thus be considered as the basis for our definitions of the spectral distortion as pure 
redistribution of photons, and in particular of the 2 : 1 division of initial energy in temperature 
perturbations into a blackbody part and a spectral distortion part. Further justification is provided 
by the fact that with this definition, the spectral distortion can also be identified with the entropy 
production. 

One important difference from the heating of CMB usually considered, for example, in clusters 
of galaxies or due to decay of particles in the early Universe, is that in mixing of blackbodies, we 
are also adding photons (compared to a blackbody at initial average temperature). The additional 
photons are able to create a new blackbody at a higher temperature, which is impossible if the 
photon number density is kept constant. 

Although we have derived the above formulae by considering only two blackbodies, they are 
applicable to an ensemble with arbitrary number of blackbodies by just replacing the average over 
two blackbodies (AT/T) 2 with the average over the ensemble , ((AT/T) 2 J. 

4. Example: Mixing of the blackbodies in the observed CMB sky 

One obvious and cleanest source of blackbodies of different temperatures is the CMB sky. The 



Komatsu et al 



tempe rature in different directions in the sky differs by ~ 10/zK dBennett et alll996 ; 
201 ll) corresponding to the 1CT 5 spatial fluctuations in the energy density of radiation/matter in 
the early Universe before recombination^. Thus any telescope, due to finite width of its beam, 
looking at the mi crowave sky will inevitabl y mix the spectra of blackbodies of different temperature 



within the beam dChluba & Sunvaevl2004i . In addition we may explicitly average the intensity over 



the whole sky to achieve higher sens itivity and precisio n in the measurement of CMB spectrum 



as is done, for exa mple, by COBE (IFixsen et al.lll996l) and in the proposed experiment PIXIE 



Kogutetal 201 1 ). The (angular) average amount of energy in CMB anisotropies is given by 



(Chluba& Sunvaev 2004) 



AE l(AT\-\ 6 v, , 

where E y is the energy density of CMB photons. One third of this energy creates a y-distortion 
of magnitude yy = j = Tilf — 8 x 10~ 10 . The measured temperature from averaged 



4 On very small scales the fluctuations are considerably smaller since these fluctuations were erased due to 
diffusion of photons before and during recombination (Silk damping). 
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intensity is also higher than the averaged temperature 7cmb by 7cmb ^(t") ^ = 4.4nK accounting 
for the remaining 2/3 of energy in anisotropics. The increase in entropy in this mixing is Scmb x 



((¥ 



3.5 x KT 6 /cm 



In the above estimate we ignored the dipole which has been measured by COBE and WMAP 



(IBennett et al 



1996; 



Jarosik et al 



20111) to have an amplitude equal to 3.355 + 0.008 mK cor- 
responding to our peculiar motion v = 1.23 x 10~ 3 . The average power from dipole is then 
3Ci/(47r) = v 2 /3. The resulting y distortion is v 2 /6 = 2.5 x 10~ 7 and increase in monopole temper- 
ature is v 2 /3 = 1.4/zK. The increase in entropy from the mixing of the CMB dipole on our sky is 
AS = 1.1 x 10- 3 /cm 3 . 



5. Application:Dissipation of sound waves in the early Universe 



Before recombination, we have a tightly coupled plasma of radiation-electrons-baryons in the early 
Universe. At high redshifts, both the energy density and pressure in the plasma are dominated by 
radiation while at low redshifts, but before recombination, the baryon energy density becomes im- 
portant, although pressure is still dominated by radiation. Sound speed in this relativistic plasma is 
therefore 1 / V3 and Jeans scale or sound horizon is particle horizon / V3. Primordial perturbation 
on scal es smaller than the Jean s scale or sound horizon theref ore oscillate setting up standing sound 
waves ( Lif shitzll 1 946 , see also Sunvaev & Zeldovichll 1970b ). Although the photon mean free path 



due to Compton scattering on electrons is very small, they are still able to traverse considerable dis- 
tance since the big bang, performing a random walk among the electrons. This diffusion and mixing 
of photons as a result of Thomson scattering erases the sound waves on scales corresponding to 
the diffusion scale (and smaller). Macroscopically, the dissipation of sound waves can be identified 
as due to the shear viscosity and thermal conduction in the relativistic fluid composed of baryons, 
electrons and p hotons. The dam ping of sound waves on small scale s due to thermal conduction was 
pointed out bv iLifshitzl (119461) and first calculated by dSilkll 19681) and is known as Silk damping. 
At high redshifts (z » 693) when the energy density of the plasma is also dominate d by radiation 
shear viscos ity is more im portant than thermal conductivity and was calculated by 



dl970l) . later lKaiserl (11983b included the effect of pola rization (see also 



ing spectral distortions in CMB were considere d by 



Sunvaev & Zeldovich 



Hu et al 



(11994a) 



berg[l< 



Peebles & Yu 



(1970a); 



Dalvl dl99ll) : 



Sunvaev & Zeldovich! dl970al) demonstrated that the upper limit to the /i-type dis- 



tortions allows us to constrain the amplitude of the primordial fluctuations, which were completely 
damped in the CMB (on small scales), and are today in the unobservable part of the matter/CMB 
power spectrum. 

Microscopically, diffusion of photons mixes photons from different phases of sound waves 
which have different temperatures. This is shown schematically in Fig. [3] This implies that locally 
a y-type distortion is created, which quickly comptonizes to a //-type distortion at z > 10 5 , as 
calculated in the previous sections. The dissipation of sound waves is best understood in Fourier 



space, denoting Fourier transform of ^-(x, n) with ©(k, n), n is the photon directioqj, x is the 



5 Bold letters denote vectors, bold letters with hat denote unit vectors and normal letters denote magnitude 
of the vector 
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X (Diffusion scale) 



Fig. 3. Cartoon picture of mixing of blackbodies: Photons from different phases of sound wave 
mix and create spectral distortion/entropy. 



comoving coordinate and k is the comoving wavenumber (see also 



Chluba et a l. 2012b), 



AE 



= 6 /(^ c """/^ <e(k ' fl,s<k - k '' fl,> 

r k 2 dk \^ 

= 6 -z-fPm ya{+\)®](k) 



(16) 



where we have expanded the temperature perturbation transfer functions in Legendre polynomial 
basis, 0(n.k, k) = (2€ + l)^(n.K)®<(A0> n is the photon direction, k is the unit vector 
along the Fourier mode which is parallel to the electron peculiar velocity in linear theory. This 
transformation is possible since in first order perturbation theory the photon transfer functions 
depend only on n.v or in Fourier space on n.k and not on n and k separately. Pj{k) is the power 
spectrum of initial curvature perturbations (£) with respect to which the transfer functions ®f 
are calculated. We have also used homogeneity and isotropy of the Universe to carry out one of 
the integrals and angular part of the second integral. Cosmological perturbations, in particular 
monopole and dipole depend on the choice of gauge. We will use conformal Newtonian gauge 
from now on. The spectral distortions defined as pure redistribution of photons, for example the y 
distortions we are considering, are however gauge independent. 

Before cosmological recombination starts with helium recombination at z ~ 6000, elec- 
trons/baryons and radiation are tightly coupled and I > 2 modes can be neglected. Most of the 
energy of sound waves is in monopole and dipole and can be calculated using relation between 
monopole and dipole in the tight coupling regime, ©i ss ©o/ V3 (using sound speed c s ~ 1/ V3), 



AE 

~eZ 



= 6 |^)K + 30|] 
C k 2 dk r , n 



(17) 



This result i s 9/4 times the estimate used in 



Sunvaev & Zeldovich 



Khatri et al 



(Il970al) : 



Hu et al 



dl994al) : 



(1201 II) where it was also assumed that all of the energy in sound waves gives rise 
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to spectral distortions. As discussed above, 1 /3 of this energy, when dissipated due to mixing of 
blackbodies, sources spectral distortions which are created as y-type but rapidly comptonize to 
type distortions or Bose-Einstein spectrum at high redshifts z > 10 5 . The remaining 2/3 of the 
dissipated energy raises the average temperature of CMB which is not directly observable. Thus 
the correct result for distortions differs from earlier estimates only by a factor of 1 /3 x 9/4 = 3/4 
Jchluba et al.lboi2bl) . 

At redshifts 10 5 < z < 2 x 10 6 , the average /i distortion therefore increases at a rate, 

ck/ _ 4 d_l AE_ 
dr ~ ~ ' dr3 ~eZ 



(18) 



and rate of increase of entropy density is given by, 



dAS/S 
dt 



1 d AE 
~4dt ~EZ 



d rk Z dk r ,1 

5 J^)[3^]. 



(19) 



The above rates are easily calculated using analytic tight coupling solutions given by 
Hu & Sugivamal (1 19951) and for power spectrum with constant scalar index it is p ossible to do 
the tim e derivatives and integral analytically. In particular the expressions presented in 



Khatri et al 



d201 lb for yu type distortions remain valid after multiplication by a factor of 3/4. The y- type dis- 
tortions at z < 10 4 require inclusion of additional modifications due to breakdown of tight coupling 
during recombination , second order Doppl er effect, and higher order temperature anisotropies, 
which were derived in 



Chluba et al 



(12012bl) using second order Boltzmann equation. 
We can use the first order Boltzmann equation, d0/df = © - ±[0 2 + ©JJ + ©^^(n.k) - 
0(n.k) - /v!Pi(n.k), to calculate the time derivative of Eq. (fT6l l. Taking into account that only 
1/3 of dissipated energy leads to spectral distortions, and requiring that the dipole/velocity term is 
gauge invariant gives us the full result]^ 



d AE 

df Ey distortion 
ignore metric perturbations 



'Ar d AT 
~dt~ 



> 4« e O"x 



Ck 2 dk 
j^ 1 



Pi{k) [©! (30! - V) 



9 -®l- l0 2 (0P + 0P) + £(2^+l)0 



impose gauge invanance 



> 4n e o"T 



k 2 dk 



2jt 2 



Pi{k) 



(30! - v) 2 



e>3 



(20) 



where \ e (k) = -i^v(k) is the transfer function of peculiar velocity of bar yons/electrons and 
denote polarization multipole moments. This expression was first derived by 



Chluba et al 



(2012b) 



We have ignored the gravitational potential/metric perturbations since they cann ot create spectral dis 



tortions. They do cancel out explicitly in the second order Boltzmann equation dChluba et al 



2012b). 



Gravitational potential/metric perturbations do contribute to the average CMB temperature but this effect 
is unobservable. 
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using the second order Boltzmann equation, which automatically takes care of the gauge indepen- 
dence and metric perturbations. Using that the spectral distortions (defined as a pure redistribution 
of photons) are gauge invariant, together with the fact that the only physical mechanism operating 
here is the mixing of blackbodies, allows us to derive the full result without referring to the sec- 
ond order Boltzmann equation and using only the well studied first order Boltzmann equation. The 
identification of a symmetry in the probl em, i.e. gauge invariance, allows us to derive very simply 
the results of the extensive calculation of Chluba et al 



(I2012bl) corresponding to the average spec- 
tral distortions in CMB created by the dissip a tion of sound waves. Using the Boltzmann equation, 
on the other hand, also allows 



Chluba et al 



(12012bl) to obtain new results on anisotropies of the 
spectral distortions, and we refer to that work for a detailed discussion. 

We note that with our definitions, all the photon transfer functions and v are real quantities. 
We have also the introduced multipole moments of degree of polarization, P , defined in the same 
way as the temperature multipole moments^ The polarization terms are coming directly from the 
first order Boltzmann equation for temperature. They contribute a t the level of ~ 5% — 10% to the 
effective heating rate close to the recombination epoch at z - 10 3 ( Chluba et alj|2012b ). 

At lower redshifts, baryons and photons develop relative velocity and second order Doppler 
effect also contributes to the y distortions and appears above in the gauge invariant combination 
with photon dipole. This effect can thus be considered as the mixing of dipole in the electron rest 
frame but in a general f rame like conformal Newtonian gauge it originates in the Compton collision 



term dHu et al 



1994bl) . This is the only significant contribution from the second order Compton 



collision term to the spectral distortions (in addition of course to the Kompaneets term). This term 
can also be easily obtained by taking the part of the second order Compton collision term (see for 
example iHu et al.lll994bt iBartolo et al.ll2007t IPitrou 2009) with the v-type spectral dependence. 



Higher order corrections originating in the terms which are second order in pertu rbation theory 
and also second order in energy transfer were calculated in IChluba et al.l (12012bl) and shown to 
be negligible. Fitting formulae for fi distortions as a function of spectral in dex and its running for 
primordial adiabatic pertu rbations are also given in 
modes were considered by ] 



Dent et al 



(2012) 



Chlub, 



Chluba et al 



from some exotic models for small-scale power spectrum. IPaier & Za ldarriaga ( 20 1 2 ) have pointed 



et al 



(l2012bl) . Recently isocurvature 



(2012a) have calculated the distortions 



Ganc & Komatsuld2012l) 



out the possibility of constraining non-gaussianity using /i distortions, and < 
have investigated the consequences of modified initial state for single field inflation. 

Equation ( f2Qb is explicitly gauge invariant and is recommended for calculations of distortions 
instead of taking the time derivative of monopole, Eq. ( fTST l. In particular, Eq. ([18} is accurate only 
in the //-era (z > 10 5 ), and the early stages of the y-era (z > 10 4 ), when only the shear viscosity 
(quadrupole) term contributes, and cannot be used to estimate the y-type distortions created at late 
times, around and after recombination. In the y-era (z < 10 4 ), thermal conductivity (dipole/velocity) 
term and I > 2 anisotropies contribute at a significant level and the full Eq. (l20b must be used. 

The first three terms in Eq. (l20l give the dominant contribution to the dissipation of sound 
waves. The first term mixes the blackbodies in the dipole resulting in transfer of heat along the 
temperature gradient, and can thus be identified as the effect of thermal conductivity. The second 
term in Eq. (l20l i. similarly, mixes the blackbodies in the quadrupole or the shear stress in the 



Our convention is same as that of IIMa & Bertschingei 



Zaldarriaga & H arari 



1995 



Dodclson 



2003) but differs from that of 



( 1995) by a factor of (—if in the definition of multipole moments. 
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photon fluid and can thus be identified as the effect of shear viscosity. The third term takes into 
account t he polarizatio n dependence of the Compton scattering and is a correction to the shear 
viscosity (lKaisenll983h . The i > 3 multipoles are negligible during tight coupling by definition 
(and thus for the //-type distortions) but give a small contributi on during recomb i nation as the 



tight co upling breaks down and the photons begin to free stream (iKhatri et al 



2011 



1995 



2012b). O n substituting the conformal Newtonian gauge tight coupling solutions dHu & Sugivama 



Chluba et al 



Zald arriaga & Hararilll995i lDodelsonll2003l) . we get for the first term (ignoring the phase of 



the oscillations which actually differs between the left hand side and the right hand side by 7r/2 

(30i - v) 2 R 2 k 2 



1 +R(n e cr T ) 2 



e?. 



(21) 



where R = jjk. — |22 anc j p b j s the energy density of baryons. At z > 10 5 , when the //-type 
distortions are created, R 2 <k 1 and thermal conductivity contributes negligibly to the sound wave 
dissip ation. The dominant terms du ring the //-type era are the second and the third (shear viscosity) 
terms ( Zaldarriaga & Hararilll995 1 



V i 02 ( p +0 p) = 



16 



(22) 



2 z 2 *\ L ' "> 15(« e o- T ) 2 
If we ignore polarization, the factor of 16/15 in the above equation would be replaced by 8/9. At 
redshifts z > 2 x 10 6 the distortions are exponentially suppressed due to the combined action of 
bremsstrahlung and double Compton, which can create and destroy photons at low frequencies, 
and comptonization, which redistributed the photons over the entire spectrum creating a Bose- 
Einstein spectrum. The rate of energy injection in Eq. (l20l therefore has to be multiplied by a 
suppression factor or blackbody visibility, Q{z), for z > 10 5 , giving the part of the energy injec- 
tio n which is actually observed a s //-type distortion. An analytic solution for Q{z) was calculated 
bv ISunvaev & Zeldovichl (11970a) . who only considered bremsstra hlung. Their solution was later 
applied to double Compton scattering bv lDanese & de Zotti ( 1982) (a ccurate to 5 - 10%) and im- 
proved recently to sub-percent accuracy in 



Khatri & Sunvaev 



(120121) . Num erical computation of 
the sp ectral distortions is possible using numerical c odes such as KYPRIX ( IProcopio & Burigana 



2009) and CosmoTher[C| (Chluba & Sunvaev 



20121) . the later code includes the energy injection 



due to Silk damping and is able to calculate such small distortions at high precision. 



6. Conclusions 

Mixing of blackbody spectra results in a photon distribution which is no longer a perfect black- 
body but contains a y-type spectral distortion. The mixed spectrum has higher entropy than the 
average entropy of initial spectra as expected from an irreversible process which creates disorder. 
The energy which goes into the spectral distortion and the increase of entropy can be calculated 
very simply using statistical physics. The increase in entropy is, at first order in small distortions, 
independent of the shape of the distortion and can be calculated using the equilibrium thermo- 
dynamics formula for small addition of heat, dS = dQ/T . The part of the energy which sources 
distortions is only 1/3 of the total energy available in temperature perturbations (with respect to 

8 This does not introduce any error in the calculation of heating of the average CMB spectrum since we 
should average the sound wave over a whole oscillation. 

9 www.chluba.de/CosmoTherm 
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a blackbody at average initial temperature) and also results in increase of entropy. The remaining 
2/3 of the energy in perturbations goes into increasing the average blackbody temperature of the 
photons and can be identified with entropy conserving part of the mixing process. We have proven 
explicitly in this paper that the comptonization of the initial spectrum with y-type distortion to the 
Bose-Einstein spectrum does not change this 2: 1 division of the dissipated energy into a blackbody 
part and a distortion part. From an observational point of view, we are just interested in the value of 
the observable fi, and this 2:1 division of dissipated energy allows us to compute the value of in a 
straightforward way. //-distortions are unique and very important because it is impossible to create 
them after z < 10 5 and thus probe the physics of the early Universe unambiguously, y-type distor- 
tions on the other hand are created throughout the later history of the Universe, and it is difficult to 
separate the contributions from the different epochs. 

We have an almost perfect blackbody in the Universe in the form of CMB. We apply our re- 
sults to mixing of blackbodies in the observed CMB sky. As a result of this mixing, the spectrum 
observed by a telescope with finite beam size would have inevitable y-type distortions. This effect 
must be taken into account in experiments aiming to measure CMB spectrum at high precision. A 
very important application of physics of the mixing blackbodies is in the early Universe. Before 
recombination, the tightly coupled baryon-electron-photon plasma is excited by primordial pertur- 
bations in energy density, resulting in standing sound waves. The photons in different phases of the 
sound wave have a blackbody spectrum with different temperature. Photon diffusion and isotropiza- 
tion of the radiation field by Thomson scattering mixes these blackbodies on scales corresponding 
to diffusion length. These spectral distortions measure the primordial spectrum on very small scales 
(with the smallest scales completely inaccessible by any other means), at comoving wavenumbers 
10 < k < 10 4 Mpc , and will thus provide a powerful new tool to constrain early Universe physics 
in the future. We have derived the energy release resulting from the damping of these sound waves, 
and the corresponding spectral distortions of the CMB, in a simple manner using the physics of 
mixing of blackbodies. These results and additio nal (but negligible) co rrections were calculated 



recently using second order perturbation theory in Oiluba et al.l d2012bl) . The results, for the very 
important fi type distortions, coincidentally are close to the estimates used in literature until now, 
with the main difference being a correction factor of 3/4. 
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